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The Anderson transition in a 3D system with symplectic 
symmetry is investigated numerically. From a one-parameter 
scaling analysis the critical exponent v of the localization 
length is extracted and estimated to be v = 1.3±0.2. The level 
statistics at the critical point are also analyzed and shown to 
be scale independent. The form of the energy level spacing 
distribution P(s) at the critical point is found to be different 
from that for the orthogonal ensemble suggesting that the 
breaking of spin rotation symmetry is relevant at the critical 
point. 
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Since the original work by Anderson |jj there has been 
considerable interest in the metal-insulator transition in 
disordered electron systems [|]J^]. Critical phenomena 
at the Anderson transition are conventionally classified 
according to the universality class into which the system 
falls: orthogonal, unitary or symplectic The Ander- 
son transition has been studied intensively numerically, 
analytically and experimentally. Nevertheless, in spite of 
this effort, we believe it is fair to say that some aspects 
of the Anderson transition remain puzzling. For exam- 
ple, the critical exponent v of the localization length has 
been estimated for three-dimensional (3D) systems with 
orthogonal || and unitary symmetry by means of 
finite-size scaling. The estimated numerical values of the 
critical exponent for these two universality classes have 
turned out to be rather close to each other. At present 
it is unclear whether or not this is an accidental coinci- 
dence. An obvious question immediately presents itself; 
does this coincidence also occur for the symplectic sym- 
metry class? And if so does this coincidence also hold for 
other characteristics of the critical point? 

Systems belonging to the symplectic universality class 
exhibit the Anderson transition even in two-dimensions 
(2D) |)-12|, while systems belonging to the unitary and 
orthogonal classes do not, in general, exhibit any Ander- 
son transition in 2D. Thus we also have the opportunity, 



when studying symplectic systems, to see how the criti- 
cal behavior at the Anderson transition depends on the 
dimensionality of the system. 

The statistical properties of energy levels P Jl^JT^ ] in 
the vicinity of the Anderson transition, and at the critical 
point in particular, have attracted considerable attention 
recently ]lq-p^]. On the metallic side of the transition 
it has been demonstrated that the energy level statistics 
can be described by random matrix theory p|, p^|JT^ ]; for 
example, the spacing s between successive energy levels 
is well described by a distribution function P(s) which 
is quite close to the Wigner surmise. In contrast, on the 
insulating side of the transition, energy levels are uncor- 
related and the spacing distribution is Poissonian. At 
the critical point, where the metal- insulator transition 
takes place, it has been claimed [ |0i||T6| ] that the energy 
level statistics are also universal but different from those 
predicted by random matrix theory. These so called criti- 
cal statistics are believed to be universal in the sense that 
they depend neither on the system-size nor on the details 
of the model Hamiltonian. It is expected that these criti- 
cal statistics should be determined only by the symmetry 
of the system and that they should be reflected in the 
critical behavior at the Anderson transition (2^,^(|. The 
existence of universal critical statistics has been demon- 
strated numerically for a 3D orthogonal system and for 
several 2D systems Jl|^-||,|3 • 

For 2D symplectic systems, it has been observed that 
the level spacing distribution at the critical point (P c (s)) 
exhibits a power-law behavior as P c (s) ~ s 4 for s -C 1 



1| |23j. For 2D unitary systems it has been found that 
the spacing distribution behaves as P c (s) ~ s 2 for s -c 1 
p2{ . These behaviors of the level correlations for small 
s are consistent with random matrix theory. However, 
it has also been reported |24|] that in 3D P c (s) is unaf- 
fected by the presence of an Aharonov-Bohm flux (AB 
flux), even though the AB flux should break time rever- 
sal symmetry. If, as is claimed in p4j , P c {s) is insensitive 
to the breaking of time reversal symmetry it immedi- 
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ately prompts the question; is P c (s) also insensitive to 
the breaking of spin rotation symmetry? 

In this letter, we consider a 3D system with symplectic 
symmetry. In particular, we estimate the critical expo- 
nent v for the localization length and examine the spac- 
ing distribution P c (s) of the energy levels in the critical 
regime. 

The model we adopt here is described by the Hamilto- 
nian |E8| 



H 



E 

<i,j> ,a,a' 



Vi,(r;j,<T'C ii(T Cj t<T i + ^ gjCj <0 .Cj )tr (1) 



with 



V[exp(-i6»cr fe )] c 



x,y,z, 



(2) 



where C\ a {Ci. a ) denotes a creation (annihilation) op- 
erator of an electron at the site i with spin a, and 
{<7 X , <J y , a z } denote the Pauli matrices. The lattice sites i 
are supposed to lie on a simple cubic lattice. The strength 
of the hopping amplitude is denoted by V. The site- 
diagonal potentials {ei} are assumed to be distributed 
independently, and their distribution is taken to be uni- 
form in the range [-W/2, W/2]. The parameter 8 charac- 
terizes the strength of the spin-orbit coupling. For 8 = 
orthogonal symmetry is recovered. Here, two cases of 
strong spin-orbit coupling, 8 = tt/6 and 8 = 7r/4, are 
studied. The length is scaled by the lattice spacing and 
x(y, z) denotes the unit vector in the £-(?/-, z-)direction. 
We note that this model is a 3D generalization of the 
Ando model JlO| proposed for a 2D symplectic system. 
It is easy to verify that the Hamiltonian is invariant un- 
der time reversal. In the following, we confine ourselves 
to the case of the band center (E = 0), for simplicity. 

First, we carry out a finite-size scaling analysis 
of the Anderson transition and estimate the critical ex- 
ponent v. We consider a quasi-one dimensional system 
[M x M x L) with L ^> M whose cross section is com- 
posed of an M x M two-dimensional lattice with a pe- 
riodic boundary condition. We calculate the localiza- 
tion length £m along such a quasi-one dimensional sys- 
tem with the transfer matrix method, in which we adopt 
the quaternion-real representation of the Hamiltonian so 
as to carry out numerical multiplications efficiently. The 
system-sizes treated here are M = 6,8, 10 and 12 with 
L up to 5 x 10 4 . An average over four independent re- 
alizations of random potentials has been performed in 
order to reduce the error of raw data, resulting in rel- 
ative errors smaller than 1%. In Fig. 1, the renormal- 
ized localization length Am defined by A m = £m/M is 
plotted as a function of the disorder parameter W for 
8 = tt/6. It is clear that the Anderson transition takes 
place at W c /V = 19.0 ± 0.2 for this 8. The standard 
analysis for the critical exponent v of the localization 
length defined as ~ \W - W Q \- V {W W c ) yields 

v = 1.3 ± 0.2. This value of v is close to those obtained 



for the orthogonal and the unitary symmetry classes in 
3D ||||-|8). The value of the renormalized localization 
length at the critical point A c is A c = 0.56 ± 0.02 (Fig. 
1), slightly smaller than those obtained for the orthog- 
onal and the unitary classes ||,[7). A similar calculation 
for 8 = tt/4 leads to W c /V = 19.9±0.2 with v = 1.3±0.2 
and A c = 0.55 ± 0.02. The fact that both v and A c are 
insensitive to the strength of the spin-orbit coupling 8 
supports the single-parameter scaling behavior of £m/M. 
By introducing the scaling parameter £, we have demon- 
strated that £m/M is indeed a single parameter function 
of M/( (Fig. 1, inset). 

Now we focus on the level statistics at the critical 
point. Clarifying the relationship between the statistics 
at the critical point and the critical behavior of the lo- 
calization length and of the conductivity is one of the 
important open problems in the theory of the Anderson 
transition. It has been proposed in refs. S| and 26 that 
the asymptotic form of the critical spacing distribution 
P c (s) of neighboring energy levels is related to the critical 
exponent v. The proposed form for P c (s) is 



P c (s) ~ exp(— Cs 



S » 1, 



(3) 



where 7 is related to the critical exponent v of the local- 
ization length 



_ ^ 1 
^ dv 



(4) 



Here C is a positive constant, which may depend on the 
symmetry of the system, and d denotes the dimension- 
ality. For s « 1, the form of P c (s) is expected to be 
P c (s) ~ s' 3 , where j3 is equal to 1, 2 and 4 for the orthog- 
onal, the unitary and the symplectic universality classes, 
respectively. It has also been suggested that the overall 
shape of P c (s) can be described by the form |pq| : 



P c (s) = As 13 exp(-Bs 2 -~<) 



(5) 



The constants A and B are determined by the two con- 
straints J P(s)ds = 1 and J sP(s)ds = 1. For a 3D sys- 
tem with orthogonal symmetry the above formula (|^), 
with v w 1.3, describes reasonably well the numerical 
data for P c (s) @|1|3- 

We have investigated the level spacing distribution 
P c (s) at the critical point for our 3D symplectic sys- 
tem, by numerically diagonalizing the Hamiltonians of 
the M x M x M systems with M = 6, 8, 10 and 12. The 
numbers of realizations of random potentials are 2000 for 
M = 6 and 1000 for M = 8, 10 and 12. One-tenth of the 
unfolded spectra around E = (\E\ ^ V) is used for cal- 
culating the spacing distribution P c (s). This is justified if 
W c does not vary considerably for energies E close to the 
band center. In fact, the critical disorder W c for E = 2V 
and 8 = n/6 is estimated to be W c /V = 18.9±0.5, which 
is almost the same as that for E = 0. 
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The spacing distributions P c (s) at the critical point 
W c = 19V (6 = tt/6) for M = 6, 8, 10 and 12 are shown 
in Figs. 2 and 3. It is seen that P c (s) is independent of 
the size of the system M confirming the existence of a 
critical spacing distribution. 

First we consider the limiting behavior of P c {s) as s — > 
0. Here we find P c (s) ~ s 13 with (3 ~ 4 (see inset Fig. 2). 
This quartic limiting behavior is characteristic of systems 
with symplectic symmetry and is a clear indication that 
the breaking of spin rotation symmetry is relevant at the 
critical point. 

Next we consider the behavior of P c {s) in the limit 
s > 1. A fit to (g) for s > 2 yields 7 = 1.0 ± 0.15. 
Making use of (Q) this yields v > 2.2 which is inconsistent 
with the value of v = 1.3 ± 0.2 we obtained earlier. In 
fact the data are well fitted by a simple exponential law 
P c {s) ~ exp(— as) with a « 1.7 (Fig. 2, solid line). We 
note that a good fit to data for 3D orthogonal Q and 2D 
symplectic systems f|l],^3| has also been obtained with a 
simple exponential. 

Finally we consider a fit to the entire spacing distribu- 
tion. If we take equation (||) with [3 — 4 and 7 = 0.75, a 
value obtained by substituting v ~ 1.3 into (Q), we find 
a very poor fit to P c (s) (dashed line Fig 2). On the other 
hand if we vary 7 so as to obtain the best fit (see Fig 
3) we find 7 = 1.43 ± 0.01 which corresponds to v < 
which is physically unacceptable. We conclude that the 
proposed form, corresponding to equations (Q) and (^) , is 
not consistent with our numerical data. For comparison, 
we also plot in Fig. 3 P c (s) for a 2D symplectic system 
[21 ^3| where the exponent v is estimated to be v sw 2.7 
[ 1 1 1 . It is worth noting that although the estimated crit- 
ical exponent v in 2D is roughly twice our estimate for 
v in 3D the corresponding values of 7 obtained from (^) 
are close to each other. If the formulae (||) and (||) were 
valid, P c (s) in 2D would have to look similar to that in 
the 3D system. We see in Fig. 3 that this is clearly not 
so. 

In summary, we have analyzed the Anderson transition 
in a 3D system with symplectic symmetry. The critical 
exponent v of the localization length is estimated to be 
v = 1.3 ± 0.2, which is rather close to the values found 
in 3D unitary and orthogonal systems. On the other 
hand, we have demonstrated that the energy level spacing 
distribution P c (s) at the critical point is sensitive to the 
breaking of spin rotation symmetry. In particular for s <C 
1 we find P c (s) ~ with f3 ~ 4 which is characteristic of 
the symplectic symmetry class. The sensitivity of P c {s), 
at small s, to the breaking of spin rotation symmetry 
is consistent with the conventional classification of the 
critical behavior according to symmetry. 
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FIG. 1. The renormalized localization length Am as a Sanc- 
tion of the disorder W for 9 — n/6. Open triangles, open di- 
amonds, open squares and open circles correspond to M = 6, 
8, 10 and 12, respectively. Inset: One-parameter scaling be- 
havior of Am for 9 — it/6 and 7r/4. Filled triangles, filled 
diamonds, filled squares and filled circles represent the data 
in the case of 9 = it/4: for M — 6, 8, 10 and 12, respectively, 
while open marks are for 9 = n/6. 



FIG. 2. The spacing distribution P c (s) at the critical point 
W = W c {= 19V) for 9 — 7r/6. Open triangles, open dia- 
monds, open squares and open circles stand for the data for 
M = 6, 8, 10 and 12, respectively. Dashed curves represent 
the result by the formula (^) with 7 = 0.75 and /3 = 4 . Solid 
line is a function oc exp( — 1.68s). Inset: Double logarithmic 
plot of P c (s) for s<l. The solid line expresses the behavior 
oc s 4 whereas the dashed lines express oc s. 



FIG. 3. The spacing distribution P c (s) at the critical point 
W — W c (= 19V) for 9 — 7r/6 is shown in the linear-scale. 
For comparison, the data of the critical spacing distribution 
for the 2D symplectic system, which is taken from Ref. p2| , 
are also plotted (crosses). The dashed curve represents the 
Wigner surmise for the Gaussian symplectic ensemble. The 
solid one represents the best fit based on the formula (|B|) with 
7 = 1.43 and (3 = 4. 
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